Landau dynamics of a grey soliton in a trapped condensate 
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It is shown that grey soliton dynamics in an one-dimensional trap can be treated as Landau 
dynamics of a quasi-particle. A soliton of arbitrary amplitude moves in the trapping potential 
without deformation of its density profile as a particle of mass 2m. The dynamics in the local 
density approximation is shown to be consistent with the perturbation theory for dark solitons. 
Dynamics of a vortex ring in a trap is discussed qualitatively. 
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Following the experimental observation of grey solitons 
(or more precisely, entities associated with quasi-one- 
dimensional grey solitons) in Bose-Einstein condensates 
(EEC's) , a great deal of attention has been paid to the 
theory of the phenomenon (see e.g. 0,13 and references 
therein). It was found that a grey soliton in a parabolic 
trap displays a number of peculiarities. Among them, we 
mention the frequency of oscillations, which is times 
less than the frequency of oscillations of the condensate 
as a whole 0, and the consequent beatings, which can 
be observed in the long-time dynamics of a soliton, an 
internal mode accompanying soliton dynamics, nontriv- 
ial phase changes in soliton evolution Q, etc. It turns 
out, however, that the mathematical treatment of the 
problem, based on the application of perturbation the- 
ory to dark solitons ^ and covering all the main effects, 
is rather involved, and more important, requires a small 
soliton velocity. The major problem is that a grey soliton 
even in the simplest one-dimensional (ID) parabolic trap 
is dramatically different from the standard grey soliton 
known from soliton theory. Distinctions emerge from dif- 
ferent boundary conditions, which are zero in the case of 
a trap potential and nonzero in the case of a mathemat- 
ical grey soliton (see e.g. 0). 

The present paper aims to describe the motion of a 
grey soliton over a wide range of velocities in a trapped 
condensate, whose longitudinal size is sufhciantly large. 
The phenomenon is described by the mean-field ID 
Gross-Pitaevskii (GP) equation 
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We consider a cigar-shaped condensate with the fol- 
lowing relations among the transverse a±, the longitu- 
dinal 2L dimensions of the condensate and a healing 
length ( = (c hereafter being the sound velocity): 

aj_ <C 1^ ^ 2L. The ID couphng constant gi, comes from 
the "averaging" over the transverse cross section of the 
condensate 0,11]. 

We developed a full description of the dynamics of the 
soliton, considering it as a quasiparticle and using the 



Landau theory of superfluidity. Subsequently we show 
that this approach is consistent with the GP equation 
and derive the main result by use of the proper pertur- 
bation procedure. In conclusion, on the basis of the re- 
sults obtained we discuss qualitatively the dynamics of a 
vortex ring. 

Soliton in a homogeneous condensate. We start with 
a gray soliton in a uniform ID condensate, i.e. when 
U{x) = 0. According to Tsuzuki the respective con- 
densate wave function can be presented in the form (see 
also i, §5.5): 



* (a:, t) = V" + - tanh 
c c 
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where X{t) — vt, v is the velocity of the soliton, n is the 
unperturbed linear density, /i — gin = mc^ is the chemi- 
cal potential, u = \/ (? — , and I = h/mu is the width 
of the soliton. The energy of the soliton is computed to 
be 
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Landau dynamics of a soliton in an inhomogeneous 
condensate. Let us now suppose that the condensate 
length, 2L, is large compared to the width of the soliton: 



(4) 



Then, for such large condensates one can use semiclassi- 
cal Landau dynamics of the soliton, where the quantity 
© plays role of the Hamiltonian of quasiparticle. We no- 
tice that condition Q is stronger than the condition of 
applicability of the Thomas-Fermi (TF) approximation 
to the condensate L ^ IS 13 ■ Requirement ensures, 
also, that the velocity of motion of the condensate as a 
whole due to the soliton oscillations is small and one can 
regard the condensate as being in the rest. 

Subject to condition Q), one can use a local density 
approximation, assuming that equation (PJ for the soliton 
energy is valid in the inhomogeneous condensate, i. e., 
that c can be changed to its local value c{X), where X 
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is the position of the center of the soliton. In the first 
approximation the soUton wave function has the same 
form ^ , where X and v are functions of time related by 
dX/dt — v{t). The sohton motion is then defined by the 
energy conservation equation 
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which obviously implies that u is also a constant. Ex- 
pressing £ in terms of it, we finally find the equation 



f 1 -<^-'' 
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which can be solved by a simple integration. 

In order to find the distribution of the sound speed, we 
approximate the density by the TF law 



TiTpix) = (m- U {x)) /gi 



(7) 



Assuming, without loss of generality, that [/ (0) = 0, we 
find 



c{X) = giTiTF I'm = Cq — U (X) /m 



(8) 



where cq = yTi/m is the sound speed at x = 0. Substi- 
tution into @ gives 

m{^^\u{X)^m{cl-u^). (9) 

This equation describes classical motion of a particle hav- 
ing mass 2to and energy m(cQ — v?) in the potential 
U {X) It is remarkable that the soliton propagates 
through the condensate without change of its density pro- 
file. Indeed, simple calculation gives for the density per- 
turbation in the vicinity of the soliton: 



onyx) = yii{x)\ — n = 
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51 cosh^ [i {x - X)] 
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This quantity does not depend of time (at a given x — X), 
because u is an integral of motion. 

Let us apply these results to a harmonic trap U {x) = 
muj^x^ /2. Notice that ujx is the frequency of oscillations 
of the center of mass of the condensate. Now the distri- 
bution of the sound velocity is {X) = c§ (l — X^/L^) 
where 2L is the condensate length defined by — 
2ix/mul and ® takes the form = ^-u^-^X^ . 

This equation describes pure harmonic oscillations with 
frequency uOg — ^xl^f^- This frequency coincides with 
the one obtained in Rcfs. 2, 3] in a different way for a 
slow soliton with w ^ cq. The amplitude of oscillations 
is equal to the coordinate X\ of the turning point: 
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The soliton oscillates between points ±Xi, where the 
soliton velocity and central density become zero. Our 
approach requires, in addition to the general condition 
that the distance between a turning point and the 
condensate boundary is larger than the soliton width: 

L — Xl ^ £. This condition can be transformed to 

1/3 

u/cq ^ (^/L) , which indicates that velocity of the 
soliton must not be too close to speed of sound. 

To conclude this section, we must make an important 
remark. Equation Q was derived for a uniform back- 
ground in the absence of an external potential. It is not 
obvious that the energy has the same form in the pres- 
ence of trapping potential. The point is that a soliton is a 
region of decreased density and, thus contains a negative 
" number of atoms" 
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One may think that to the energy £ one must add the po- 
tential energy NsU{X) that corresponds to these atoms. 
Such conclusion would be wrong. Actually, this decrease 
of energy is compensated by increasing the energy of 
atoms outside of the soliton. We will prove this in the 
next paragraph by direct calculation of the energy. 

It is worthwhile to notice that the depletion (|12|l of 
the number of atoms in the soliton depends only on its 
energy £. This implies that Ng is also an integral of 
motion. A soliton is a quasiparticle of constant (negative) 
mass niNs- However, this "physical" mass of the soliton 
depends on energy and is not equivalent to the " effective 
mass" 2m. 

Energy of a soliton in an external field. In this para- 
graph we will show that equation (PJ for the energy of a 
soliton, obtained for an uniform condensate, is actually 
valid also for a trapped condensate in the local density 
approximation. Thus the trapping potential does not en- 
ter explicitly in the expression for the energy of a soliton. 

To take properly into account the expelling of atoms 
from the soliton, it is more convenient to work not at 
constant number of atoms, but at constant chemical po- 
tential fi. Thus, instead of the energy i?, we will consider 
the "grand canonical energy" E' = E — fi J n{x)dx. In 
the GP approximation this energy is £" = J e' dx, where 

e'(a^) = 1^1*=^!' + \9Mxf + {U{x) - ^i)n{x) . (13) 

Let the soliton center be at x = X. We separate 
the integration on two domains E' =^ S\x-x\>s^' '^^ ~'~ 
i\x-x\<& '^li^re L ^ 6 ^ i. In the first integral one 
can use the TF approximation for the energy, neglect- 
ing the gradient term: e'rpp{x) = ^ginTpix)'^ + {U{x) — 
IJ,)nTF{x), while, taking into account the smooth behav- 
ior of the potential, in the second term one can approxi- 
mate U{x) « U{X). Next, one can add to the first term 
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an integral f\^_x\<:5 ^'tf'^^ correspondingly, deduct 
it from the second term in E' . Such an addition trans- 
forms the first term in the total energy of the condensate 
in the absence of the soliton, which will be designated 



E'q. In the second term one can safely change the smooth 
function e'rpp{x) to e'rpp[X). Finally we arrive at the re- 
sult 
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+ 7;9iH^f - riTFiXf) + iU{X) ~ ii){n{x) - nTF{X)) 



dx 



(14) 



Now we can use the theorem about small increments 
(see §16). According to this theorem small correc- 
tions to E and E' , are equal, being expressed, corre- 
spondingly, in terms of N and /x. Thus to obtain a cor- 
rection to the energy, we must express the second term 
in (|14|l in terms of the density. One may eliminate ji 
from (|14|l by taking into account that, according to Q, 
(/i — U{X)) — giUTFiX). Owing to the fast convergence 
of the integral, it is also possible to change formally the 
integration limits from X ± (5 to to ±oo respectivily. Fi- 
nally we get E = Eq + £ where 



^\^i>^\^ + lg,(n{x)-nTF{X))' 
2m 2 



dx . (15) 



Equation (|15|l does not contain the trapping potential 
explicitly. Substituting '^(x) from ^ and integrating, 
we obtain equation (O, where n is changed to nTF{X). 
Equation ensures now conservation of E. 

Perturbative approach to grey-soliton dynamics. In 
the previous paragraphs investigating the soliton dynam- 
ics we postulated, in the spirit of the Landau theory, that 
when the soliton moves in a weakly-inhomogeneous back- 
ground its local energy stays constant. Let us now inves- 
tigate the relation between the above phenomenological 
approach and the evolution of a grey soliton emerging 
from the dynamical approach, based on the mean-field 
GP equation (P). 

A grey soliton propagates against a background, and 
thus the first step to obtain the soliton dynamics is to de- 
termine the background (which above was approximated 
by the TF density distribution). One can do this by 
means of the ansatz = ex.p{-i!^t)F{x)ip{x,t) where 
F{x) is a real-value solution of the nonlinear eigenvalue 
problem 



liF = - — F,, + Uix)F + gmF^ 
Zm 



(16) 



with F(0) = 1,^;(0) = and \im^-,±oo F{x) = 0. As 
above we stipulate that U{0) = and without loss of 
generality impose the normalization condition on F{x) at 
X = assuming that the soliton is placed in the interval 
X e {—Xi,Xi), where Xi ^ L. The function ip{x,t) 



then solves the equation 
Zm 



gM'-n)^j = R[^b,F] (17) 



which is subject 

lim2;^±oo 1-0(2^,^)1 = ■ 



to the boundary conditions 
/n. R[ip,F] is given by 



^ 2m F ^ 



+ 5in(F2-l)(|^|2-n)0. (18) 



An explicit form of Eg. ljlTII appears to be useful if 
F] is small compared with the left hand side of 
that equation. In order to estimate it we recall condi- 
tion , which can be viewed as a definition of the small 
parameter of the problem e — IjL. We restrict our dis- 
cussion to a soliton near the center of the trap, where 
g\n ~ 0{e) and the potential is smooth enough, 
which enables one to neglect third order in the expansion 
around a; = 0. One can ensure then that \i ~ g^n and 
i?[V', F\ = 0{egin) in the domain specified above. 

This naturally leads one to use perturbation theory, 
the first step of which would be an adiabatic approxima- 
tion, i.e. an approximation in which the soliton shape 
is given by Q with slowly varying parameters. This 
straightforward approach, however meets difficulties in 
the case at hand because of the first term in the right 
hand side of Eq. I|18|). Mathematical reasons for this 
are discussed in 0,3^ Here we will use a more physi- 
cal approach. Namely, we will prove that it is possible 
to define an adiabatic approximation for a dark soliton, 
such that in the leading order the number of particles 
associated to the soliton is constant. In other words, one 
has to prove that there exist a function (f){x,t), satis- 
fying ^ in the leading order and preserving the quan- 
tity A^s[0] — J\<j>\'^dx. Because both Ns and £ are ex- 
pressed in terms of the same quantity u, the energy 
£ is also an integral of motion in accordance with the 
Landau theory. By direct algebra, one may verify that 
the first term in the r.h.s. of (|18|l results in change of 
the number of particles in the case of a nonzero current 
[cx f '^{ipx'ip* ~ 4'i'x)dx]. In order to apply our argu- 
ments, we again assume that in the adiabatic approxi- 
mation the soliton is given by Q. One can then ensure 
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that the imaginary part of all terms in H17() becomes ore- 
cisely zero in the vicinity of the point x = X{t) except 
for the first term in the r.h.s., which even grows with 
X{t). This behavior originates from non-adiabatic ef- 
fects, which must be taken into account before one passes 
to the Landau quasiparticle description the soliton. 

In order to avoid this difficulty, let us take into account 
that the solution we are interested in is a function of t 
and C ~ X — X{t), where X{t) is the coordinate of the 
soliton center and can be associated with a coordinate in 
the Lagrangian description of the condensate flow. Then, 
we make the substitution 

i;{x,t)=^{x,t)-if{t)cj,^{x,t) , (19) 



2m Zm 



Now one can verify that dNs[(l)]/dt — 0{e^), i.e 4>{x,t) 
describes a soliton with a constant mass. Indeed, while 
the last two terms in the r.h.s. of 121|) obviously give zero 
contribution to the change of the number of particles, in 
order to treat the first one we recall the above require- 
ment of the smoothness of the potential U{x). One then 
easily estimates 

FA^) _ F^Xjt')) ^ {x~X)U,.M , 2^ 
F{x) F{X{t')) 2gn ^ ^ 

It is significant that the r.h.s. is an odd function of a; — 
X. On the other hand (t>*(f>x ~ 0x'/'j with (/) given by the 
r.h.s. of ©J is an even function. Thus, in computing 
dNs[(l)\/dt, the first term contains an odd function of a: — 
X , which yields zero to the leading order. 

According to the previous reasoning, the result ob- 
tained proves the conservation of the soliton energy. 

Vortex ring. The method developed in this letter can 
be applied to other localized excitations. The most inter- 
esting example is a vortex ring. Let a small ring moves 
along the axis of the trapped condensate. The energy 
of the ring in an uniform condensate can be written as 
£ = y/nx{v/c), where function x is a decreasing function 
of its argument. (It was calculated in see also @, 
§ 5.4.) Then it follows from conservation of £ that when 
the ring moves in the direction of decreased density, the 
ratio v/c (and v itself) decreases and its radius increases. 
On the contrary, moving in the direction of increasing n, 
the ring accelerates. If the central density is big enough, 
the ring velocity can reach its maximal value v = 0.93c. 
The behavior of the system near this point cannot be de- 
scribed in the local density approximation. Probably the 
excitation collapses beyond this point. 



where 



m)-l r^^WOldt' (20) 

■^^''-mJo F{X{t')) ■ ^'"^ 

Here (p can be associated with the soliton wave function 
while the second term on the r.h.s. of Eq. (|19() is an in- 
ternal mode excited when a soliton moves in a potential. 
Neglecting the terms of orders higher than e one obtains 
the equation for 0(x, t) 



<t>x + gMF'' - l)m^ - n)<l, - 2if{t)g,^^cf>l . (21) 

I 

In conclusion, we have shown that considering a grey 
soliton as a quasiparticle in the spirit of the Landau the- 
ory of superfluidity, one can obtain a simple solution of 
the problem of soliton motion in a trapped ID conden- 
sate. The energy and the shape of the soliton are pre- 
served during its motion and soliton moves in a trapping 
potential as a particle of mass 2m. 
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